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Temporal fluctuations of the atmospheric piston are critical for interferome- 
ters as they determine their sensitivity. We characterize an instrumental set- 
up, termed the piston scope, that aims at measuring the atmospheric time 
constant, r , through the image motion in the focal plane of a Fizeau interfer- 
ometer. 

High-resolution piston scope measurements have been obtained at the observa- 
tory of Paranal, Chile, in April 2006. The derived atmospheric parameters are 
shown to be consistent with data from the astronomical site monitor, provided 
that the atmospheric turbulence is displaced along a single direction. 
Piston scope measurements, of lower temporal and spatial resolution, were for 
the first time recorded in February 2005 at the Antarctic site of Dome C. Their 
re-analysis in terms of the new data calibration sharpens the conclusions of a 
first qualitative examination [1]. 
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1. Introduction 

Interferometers have been introduced in astronomy to gain in spatial resolution without the 
need to build extremely large telescopes. To resolve Sirius, observations in the infrared do- 
main (~ 2 /im) would require a telescope of about 170 m mirror diameter. Fortunately, Sirius 
can also be resolved by two telescopes of more modest size, separated by 170 m and oper- 
ated as an interferometer. Yet, despite this considerable gain in resolution, interferometers 
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are not the prime tool of today's astronomers. This is largely due to their limited sensitiv- 
ity: atmospheric turbulence makes the interferometric fringe pattern move in the detector 
plane. Accordingly, one tends to use exposure times that are short enough to "freeze" the 
turbulence, i.e. typically several milliseconds. To increase the sensitivity, phasing devices are 
being designed that measure the position of the fringe pattern due to a reference star, and 
correct continuously for the fringe motion of the target object. For such devices to work, a 
sufficient number of photons need to be collected on the reference star during the time when 
the atmosphere is frozen, i.e. during the atmospheric coherence time Tq = 0.314 Tq/Vq/s, 
where r is the Fried parameter and V5/3 is a weighted average of the turbulent layers' ve- 
locities. Clearly, the coherence time is the parameter that determines the performance of 
today's interferometers. Different definitions of the atmospheric coherence time have been 
introduced in relation to various observational techniques: single telescopes with or without 
adaptive-optics, interferometers with or without fringe trackers etc. However the standard 
adaptive-optics coherence time To has been shown to quantify the performance of all these 
techniques [2]. 

In a previous article [1], we characterized the temporal evolution of fringe motion at Dome C, 
a summit on the antarctic continent, and a potential site for a future interferometer, using 
the motion of the fringe pattern formed in the focal plane of a Fizeau interferometer. The 
temporal and spatial sampling of the measurements were low due to the available equipment 
and, instead of determining coherence-time values, the mean duration of correlation was 
assessed by fitting the fringe correlation-function onto an exponential curve (cf. Section 4). 
Such measurements have now been repeated at the site of Paranal, Chile, with sufficient 
spatial and temporal sampling, to allow the determination of the coherence time. Further, 
all relevant atmospheric parameters are constantly monitored at Paranal by a meteorological 
station, hence the parameter values derived through our set-up (termed piston scope) can 
be checked against reference values. 

In the first Section, the quantities measured with the piston scope are related to the fol- 
lowing atmospheric parameters: the Fried parameter, the turbulent layers velocities and the 
coherence time - using the Kolmogorov theory of atmospheric turbulence. The relations are 
then tested on the observations performed at Paranal. It is shown that when the sampling is 
sufficient, the precision on the coherence time is limited by the piston scope's sensitivity to 
wind direction. Given these results, the third Section presents a new analysis of the measure- 
ments obtained at DomeC [1]. The lower limits to the coherence time, derived through our 
first qualitative analysis, are confirmed and additional results on the Fried parameter and 
wavefront speed are given. 
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2. Formalism 



The purpose of the piston scope experiment is to track the rapid fluctuations of the atmo- 
spheric piston. To this effect, the entrance pupil of a telescope is covered by a mask with 
two circular openings. The resulting image is a fringe pattern within the superposition of 
the two diffraction discs. Atmospheric turbulence keeps the image of the star moving on the 
detector. The local inclination of the wave front over each of the holes causes the movement 
of the Airy discs, whereas difference in the optical path for the two holes, i.e. the piston, 
shifts the fringe pattern relative to the center of the Airy discs. Telescope vibrations, on 
the other hand, cause merely a common movement of Airy discs and fringes. The relative 
movements between Airy discs and fringes are, therefore, solely due to the atmospheric tur- 
bulence. The subsequent analysis deals with their temporal patterns. Piston changes shift 
the fringe pattern relative to the Airy discs along the interferometric axis. Accordingly, in 
order to assess the temporal fluctuations of the piston it is sufficient to consider the shift 
along the axis. 

As suggested by Conan et al. [3], the spatial power spectrum of the relative move- 
ments between Airy discs and fringes is derived from the phase spectrum assuming a 
Kolmogorov model of turbulence with an infinite outer scale. In the following we use the 
notations of Conan et al. [3]. 

/•+oo 

WJf) = 0.00969 k 2 / f~ u/3 C 2 dh, (1) 
Jo 

where f is the spatial frequency and k = 2ixj\ the wavenumber. The turbulence intensity of 
a layer % of thickness dh at altitude h is specified in terms of C 2 dh. The explicit dependence 
of C n and all following parameters on h is dropped to ease the reading of the formulae. The 
measured quantity is the separation - along the interferometric axis, x - between the central 
fringe and the center of the combined Airy discs. The spatial filter M that converts W v into 
the power spectrum equals: 

M(f) = A/(2tt) A(f)FT[(6 B -5 )/B-(5 B + 6 )/2*d/dx](f), (2) 

for a baseline vector B and the aperture filter function A(f). For a circular aperture of 
diameter D, A(f) = 2J\{'K f D) / {tt f D) and / = |f|. J n stands for the Bessel function of order 
n. FT represents the Fourier transform, 5^ is the delta function centered on L and * denotes 
convolution. Hence, 

M(f) = A/(2vr) A(f) [2 sin(TrfB) - 27rfB cos(TrfB)] / B (3) 
W^f) = M 2 (f)Wyf). (4) 

In the single layer approximation, we assume the turbulent layer to be transported with a 
velocity V directed at an angle a with respect to the baseline. The temporal power spec- 
trum of the measured quantity is obtained by integrating in the frequency plane over a line 
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displaced by / x = u/V from the coordinate origin and inclined at angle a. Let f y be the 
integration variable along this line and f 2 = f 2 + f 2 . The temporal power spectrum equals: 

1 r+oo 

w 4>( u ) = 77 W <f> (/xcosa + f y sin a, f y cos a - / x sina) d/ y (5) 

V J — oo 

We then derive the expression of the structure function: 

/+oo 
(1 - cos(27ri/t))tity(i/)dz/ (6) 
-oo 

/■+oo 

= 2 x 0.00969 C 2 dh / B 2 /~ 8/3 (2 MnfcftKTcfdWdf 

Jo 

I *(l - cos(2tt/ cos(# + a) Vt)) 

JO 

[2 sin(7r5/ cos 0) - 2nfB cos cos(tt /£ cos #)] 2 dfl (7) 

The best estimate of the parameters is obtained by fitting the measured points to: D^fy + K, 
where K is a constant that allows for white measurement noise. As seen from Eq. 7, the 
structure function depends on the wind orientation a because the mask of the piston scope 
is not rotationally symmetric. Temporal evolutions of the structure functions, for different 
values of a, are represented on Fig. 1. The asymptotic value of the structure function at large 
time increments is determined by the Fried parameter r , whereas the time needed to reach 
the asymptotic value is a function of the velocity V . 

3. Measurements at Paranal 

3. A. Observational set-up 

Several observations of Spica were obtained at Paranal on the nights from 22-23 and 23-24 
April 2006, using a modified SLODAR [4] (Slope detection and ranging). This SLODAR is 
designed to measure profiles of the atmospheric turbulence with a telescope that has a 0.4 m 
diameter primary mirror, and a focal length of 4.064 m. The detector is a 128 x 128 array of 
(24 x 24)/xm 2 pixels with a peak quantum efficiency of 92% at Ao = 550 nm and next to zero 
read-out noise. For our experiment the entrance pupil of SLODAR was covered by a mask 
with two circular openings of diameter D = 0.115 m and centers B = 0.260 m apart. The 
resulting image is a fringe pattern of angular period Xq/B = 0.44" within the superposition 
of two Airy discs of diameter 2.44Ao/-D = 2.41". Two lenses were used to increase the focal 
length by a factor 16.67, this makes each pixel correspond to an angular increment of 0.073". 
During the first night, a sequence of 1000 images was recorded at 240 Hz with an exposure 
time equal to 2 ms. On the following night, six sequences of 1000 images were recorded at 
300 Hz with 1 ms exposure time. 

The piston is quantified in terms of the motion of the fringe packet relative to the combined 
Airy discs. The quantification of the axial motion requires the extraction of the following 



4 



parameters from the observed images: the position of the central fringe and the position - 
along the interferometric axis - of the center of the combined Airy discs. This extraction has 
been described in detail in a previous article [1]. An example of a raw image is shown on 
Fig. 2 with the corresponding, fitted intensity profile. 

3.B. Derivation of atmospheric parameters 

The Fried parameter tq, the wavefront velocity V and orientation a are derived by fitting 
D^it) + K onto the data points, as described in Section 2. D${b) corresponds to an atmo- 
spheric model where the turbulence is contained in a single layer, that is displaced as a whole 
with the velocity V under an angle a. 

The resulting parameter values and uncertainties are indicated on Fig. 3. The latter corre- 
spond to a doubling of the squared deviation of the data points to the theoretic structure 
function. The Fried parameter is determined by the asymptotic value of the structure func- 
tion at large time increments. To ease the comparison with the meteorological station of 
Paranal, we indicate the seeing angle eo rather than the Fried parameter ro, these two pa- 
rameters are essentially equivalent: e = 0.976 A/r [rad]. V and a are derived from the first 
few measurement points and the coherence time, To, is then obtained through the classic 
relation: r = 0.314 r§jV . 

3. C. Performance of the piston scope 

On Figs. 4-6, the values of e , Vp S and r obtained with the piston scope are compared to 
measurements in terms of the Paranal monitoring-instruments. We do not compare the 
wind orientations, because the value of a that is obtained with the piston scope depends on 
the position of the mask, hence on the pointing of the telescope, and it is difficult to relate 
it to the angle measured by the meteorological station. 

• Seeing values (see Fig. 4): The values estimated with the piston scope coincide with 
those measured at 6m height by the DIMM [5] (Differential Image Motion Monitor). 
Note that we assume the atmosphere to consist of one layer displaced along a single 
direction, yet Fig. 7 shows that on 22-23 April the turbulence was contained in several 
layers with similar intensity. However, the asymptotic value of the structure function 
has the same altitude dependence as the seeing: 



therefore the seeing estimated by the piston scope is correct independently of turbulence 
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profile. 
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Velocities (see Fig. 5): The wavefront velocity V ps derived with the piston scope is a 
turbulence- weighted average of the layers' velocities V(h). Ideally, Vp S should have the 
same dependence on turbulence parameters as r , hence: 

V PS oc V 5/3 - [ jrc2n(h)dh ] (9) 

Sarazin & Tokovinin [6] give an empirical relation between V5/3 and the wind speeds 
measured at ground level and at 200 mB pressure. That relation has been verified 
at Paranal and Cerro Pachon in Chile, and later confirmed at San Pedro de Martir, 
Mexico: 

F 5/3 « max(^ g , 0.4 V 200mB ) (10) 

At Paranal, V g is measured by wind sensors at 30 m height and V^oomB is estimated every 
6 hours by the ECMWF [7] (European Center for Medium Range Weather Forecast) 
through a global meteorological model which runs twice a day at 00 UT and 12 UT. This 
involves the assimilation of worldwide-collected data from radio soundings, satellite 
observations etc. 

It appears from Fig. 5 that the wavefront velocities derived with the piston scope 
coincide with 0.4 V^oomB, rather than V5/3 ~ max(Vg, 0.4 V^oomB)- When the turbulence 
is contained in several layers, the measured structure function is an average of single- 
layer structure functions as represented on Fig. 1. If these layers have different wind 
velocities and orientations, the dispersion of the data points around the best-fitting 
structure function is large and the resulting wavefront velocity is poorly constrained. 
Accordingly, and in line with Fig. 7, V ps is derived with respectively 55% and 10% 
uncertainties during the first and second night of observations. 

The coherence time (see Fig. 6) is a combination of the seeing and wavefront velocity, 
thus it is essentially unconstrained during the first night. On the subsequent night, 
the values are consistent with those derived through the two following methods: With 
MASS, t is assessed from the scintillation through a 2 cm diameter aperture. MASS 
is not sensitive to the lower layers of turbulence (< 500 m), and, correspondingly, 
measures higher coherence times. A second value of To is obtained by combining DIMM 
seeing-values with measurements of the wind speed: To = 0.314 ro/^5/3, where V5/3 
is estimated by Eq. 10. Since these values are obtained from distinct locations with 
different telescopes pointing at different stars, we do not expect them to coincide. The 
results seem to suggest that the piston scope sees more turbulence than MASS and 
DIMM: While this is probable - the piston scope is installed inside a dome at ground 
level, whereas MASS and DIMM are placed on an open platform at 6 m above the 
ground - no definite conclusion is possible given the amount of data. 
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4. Measurements at Dome C 



Dome C is one of the summits on the Antarctic plateau with altitude 3233 m. The station, 
which is jointly operated by France and Italy, is located 1100 km inland from the French re- 
search station Dumont Durville and 1200 km inland from the Italian Zuchelli station. Dome C 
is known as a site with extremely low wind speeds at high altitudes. Because of its ele- 
vated location and its relative distance from the edges of the Antarctic Plateau, DomeC 
does not experience the katabatic winds characteristic of the coastal regions of Antarctica. 
Hence the coherence times could be particularly high. Lawrence et al. [8] have, during the 
Antarctic night, determined high-altitude turbulence parameters that are 2 to 3 times better 
than at mid-latitude sites. Accordingly, they concluded that an interferometer located on 
Dome C might allow projects that would otherwise require instruments in space. The value 
of To = 7.9 ms obtained by Lawrence et al. was derived from measurements with the MASS 
instrument and hence, it does not take into account the turbulence below 500 m (see the 
MASS website for corresponding calibration studies [9]). Measurements of To integrated over 
the whole atmosphere still need to be obtained. 

In this context, similar measurements to those presented in Section 3, have been performed 
at Dome C, Antarctica, on January 31 st and February 1 st 2005 at daytime. For these measure- 
ments, Canopus was observed with a telescope of focal length 2.80 m and a primary mirror 
of 0.28 m, placed 3.5 m above the ground. The entrance pupil was covered by a mask with 
two 0.06 m diameter circular openings and centers 0.20 m apart. The observational set-up, 
as well as a first qualitative data analysis has been presented in a previous article [1]. The 
observations - done with the available equipment - were both spatially and temporally under 
sampled. During six sequences out of nine, it was nevertheless possible to place a lower limit 
equal to 10 ms to the mean duration of correlation t c of the fringe patterns. This was done 
by fitting an exponential curve onto the measured structure functions: 

D^t) = D^t > t c ) x (1 - exp(-(t/t c )- 5 / 3 )) (11) 

In Section 2, the structure function has been related to the Fried parameter ro and to the 
velocity vector V in the case of a single turbulent layer, using the Kolmogorov model of atmo- 
spheric turbulence. This relation has been tested on well sampled piston scope measurements 
recorded at Paranal (see Section 3), and is now applied to re-analyze the data from Dome C. 
We consider six out of nine sequences that were presented in the previous article. Images 
were taken every 28 ms, with exposure times of 1, 2 or 3 ms. Each sequence contains between 
209 and 723 images and, thus, lasts roughly 5 to 20 s. Two sequences - recorded on February 
1 st at 7:49 UT and 9:41 UT - are not re-analyzed because the central positions of the fringe 
pattern and of the combined Airy discs are determined with too large uncertainties. In the 
previous article they were part of the three sequences during which the correlation time t c 
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was found to be less than 10 ms. For the third such sequence, recorded on January 31 st at 
9:07 UT, the fringe pattern can be fitted but since the structure function reaches its asymp- 
totic value at the first measurement point, it can not be compared to a theoretical curve. 
As seen on Fig. 8, the structure functions reach their asymptotic value after the 4 th to 5 th 
data point. The fit involves three free parameters eo, V,a besides the white noise, K, that 
is approximately constant if the instrumental settings do not vary. The data obtained at 
Paranal from April 23 rd to 24 th yield: K = (3.2 ± 0.7) 10~ 14 rad 2 . To constrain the fit, K 
is therefore fixed to the value that optimizes the global result of the six fitting procedures: 
K = 1.1 x 10~ 12 rad 2 . 

The derived parameter-values and uncertainties are indicated on Fig. 8. As specified in Sec- 
tion 3, the uncertainties correspond to a two-fold increase in the squared deviation of the 
data points to the theoretic structure function. The values of the seeing are consistent with 
measurements by DIMM (Fig. 9): The difference in the estimates by the piston scope and 
the DIMM at 8.5 m height, resembles the scatter between the values estimated by the DIMM 
instruments at 3.5 m and 8.5 m, and is due to ground layer turbulence. In line with our pre- 
vious qualitative analysis, coherence times are found to lie above 10 ms during the periods 
when five of the nine sequences were recorded. 

Note that the wind orientations are not constrained by the analysis. To derive - without 
continuous assessment of wind-direction profiles - more accurate values of tq, a parameter 
needs to be measured that is independent of the wind orientation. We have pointed out what 
appears to be a suitable new method in a previous article [2]. 

5. Conclusions 

The atmospheric coherence time, To, is the crucial parameter for interferometers because it 
determines their sensitivity. Yet, a simple method is still lacking to monitor the coherence 
time at different sites, and to decide where the future large interferometers ought to be built. 
Does the piston scope fulfill this need? To answer that question, we have related the meas- 
ured quantity to parameters of the Kolmogorov model of turbulence. 

It was found that due to its sensitivity to the wind direction the piston scope can be used to 
assess the wavefront vrelocity and the coherence time if, and only if, the whole turbulence is 
displaced along a single direction. Since the single layer model is not a permanent feature on 
most sites, the estimation of the coherence time is insecure. This conclusion is supported by 
seven sequences of 1000 images, recorded with the piston scope at the observatory of Paranal 
in April 2005. To determine the coherence time for any kind of atmospheric turbulence, a 
rotationally symmetric set-up has been proposed [2] and first measurements are planned. 
The measurements performed at DomeC have been analyzed using the method here pre- 
sented. Within the uncertainties due to low samplings, seeing angles are derived that coincide 
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with simultaneous DIMM measurements. Mean wavefront speeds are found to be remarkably 
low. In agreement with a first qualitative analysis [1], the corresponding coherence times are 
determined to be superior to 10 ms during five out of nine sequences. 
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Fig. 1. Structure functions of the fringe position relative to the combined Airy 
discs, for an interferometer with mirror diameters D and baseline length B. 
The atmosphere is assumed to consist of a single layer displaced with wind 
speed V at an angle a from the baseline. The values of a are indicated in the 
bottom right box. 





Fig. 2. Example of an image recorded with 1 ms exposure time at Paranal on 
the night of 23-24 April at 02:03:55 UT and fitted intensity profile along the 
axial direction. 
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Fig. 3. Theoretical structure functions (dashed line) fitted onto data obtained 
at Paranal, the resulting seeing e , velocity V, wind orientation a and coher- 
ence time To are indicated. 
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Fig. 4. Seeing values measured at Paranal with the DIMM and the piston 
scope. The uncertainties of the piston scope values correspond to a twofold 
increase in the quality of the data adjustment. 
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Fig. 5. Wavefront velocities obtained with the piston scope (Vp S ), wind ve- 
locities measured by sensors at 30 m above the ground of Paranal (V g ) and 
interpolated at 200 mB from ECMWF data (VboomEs)- 




Fig. 6. Coherence times obtained at Paranal through three different methods. 
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Fig. 7. Profiles of the free atmosphere turbulence obtained by MASS at 
Paranal. On 22-23 April (left panel) the turbulence was contained in several 
layers of similar intensity, while on 23-34 April (right panel) one layer at 4 km 
was predominant around 2:00 UT. 



14 



31 January 05 - 9:03 UT 



1 February 05 - 7:45 UT 



2.0*10~ IJ - 



€ = 0.52 + -0.03 " 
V= 1.73+-0-71 m/s 

a= 0-90 * 

t = 35.4 + -16.6 ms 



5.0x10" 13 - 



£ = 0.58+-0.02 " 
V= 1.98+-0.39 m/s 
a= 0-90 
t = 27.7+- 6.3 ms 



1 February 05 - 8:1 1 UT 



e D = 0.36+-0.01 " 
V= 1.88+-0.43 m/s 

a- 0-90 ° 

T n = 47,0+ -12.4 ms 



1 February 05 - 8:56 UT 









£ = 1.30 + -0.04 " 








V= 2.01+-0.59 m/s 








a= 0-90 ° 








t = 12.1 -I — 3.9 ms 













0.1 


0.2 


0.3 



1 February 05 - 8:59 UT 



e = 1.24+-0.02 " = 

V= 2.32 + -0.32 m/s - 

a= 0-90 :. 

t = 1 1 .0 + - 1 .7 ms = 



1 February 05 - 9:0D UT 



e = 1.23 + -0.02 " 
V= 2.99 + -1.28 m/s 
a= 0-90 
T = 8.6+- 3.8 ms 



Fig. 8. Atmospheric parameter values derived from measurements at Dome C. 
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Fig. 9. Seeing values measured at Dome C with the DIMM and the piston 
scope. 
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